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Abstract. Global visibility of naked singularities is analyzed here for a class of 
spherically symmetric spacetimes, extending previous studies - limited to inhomo- 
geneous dust cloud collapse - to more physical valid situations in which pressures 
are non- vanishing. Existence of nonradial geodesies escaping from the singularity is 
. shown, and the observability of the singularity from far-away observers is discussed. 
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1. Introduction 



The study of gravitational collapse of spherically symmetric solutions in General 
Relativity led to many examples of locally naked singularities, starting from pioneer- 
ing works in the early 80's - a quite exhaustive and updated list of references can be 
■ found in reviews lIHllini- Most of these papers concentrate on the aim to find photons 
CN ! emanated from the singularity and escaping the Schwarzschildian trapped region at 
^ I least locally. Moreover, the analysis has been always limited to the study of radial 
O ' null geodesies, quite simplifying then the system to study, which remains not defined 
^ ■ at the singular point, but fully decouples. 

--^ ! However, in order to study the effects of these photons for a distant observer, 

qh| global behavior of null geodesies must be studied in full generality, since nonradial 
^ ■ geodesies, with nonzero angular momentum, determine the angular diameter of the 
bX). central naked singularity as seen by the observers, giving the measure of the "size" 
of the singularity. In this direction there are in literature quite detailed studies of 
k>( ; Tolman-Bondi self similar dust cloud collapse. In [I2ll necessary conditions - 
5^ ' which actually turn out to be sufficient - for nonradial geodesies existence are de- 
rived, and behavior at infinity of these photons is numerically studied for some par- 
ticular case. Further, a complete result of nonradial geodesic existence but under the 
assumption of self similarity of the dust solution is fully proved in [|14J . and consid- 
erations about topology of the singularity are derived. 

Of course, the dust model cannot be considered as a good physical model of sin- 
gularity formation, because pressures are expected to occur during the collapse. In 
this paper we extend the results on dust models to a wide class of solutions, found 
in fT\, for which local nakedness results from a suitable choice of initial data. This 
class represents the wider set found so far of locally naked singularities in the grav- 
itational collapse of elastic materials. Since these models are not pressureless, and 
in particular radial pressure does not generally vanish along a timelike hypersurface, 
a junction between a solution from this class and a anisotropic generalization of de 
Sitter space time discussed in [SJ is performed to construct a global model. It is found 
that properties satisfied by dust solution and showed in previously cited works remains 
qualitatively valid for this wider class of collapsing spacetimes. 

The paper is organized as follows. Interior and exterior solutions are described 
in Sectional and local existence of nonradial geodesic - which are continued in the 
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exterior spacetime using junction conditions - is derived in Section 13 The resulting 
behavior for a distant observer is analyzed in Section HI together with some remarks 
about photons with infinite redshift. Last Sectionals devoted to overall conclusions. 

2. Collapse of spherical anisotropic matter 

In this part we briefly review the model we will deal with for the rest of the paper. 
The spacetime will be made by two different solutions to Einstein field equations: 
an interior part, that collapses until singularity forms, and an exterior part, that is 
matched with the former at a timelike hypersurface in order to satisfy Israel-Darmois 
conditions 1 9 1 ; it amounts to say that both first and second fundamental forms induced 
by the two solutions on this hypersurface respectively coincide. The interior part is 
described by the collapsing anisotropic solutions found in [7|, which will be matched 
with the de Sitter generalizations discussed in |5| (see also references therein). 

2.1. Interior region. The interior region is provided by a class of anisotropic col- 
lapsing matter. To describe it, the use of area-radius coordinates, first introduced by 
Ori lilSJ . turns out to be more useful than usual comoving reference frame 

(2.1) ds^ = -e^^dt^ + -dr^ + R\d9^ + sin^^df/?^). 

V 

Indeed, as well known flT^, the general spherical matter distribution can be given 
assigning energy density e and pressures Pr and pt as suitable functions of the triple 
(r, R, rj). Under the assumption ^ = 0, in [VJ it is shown that the metric takes the 
form 

(2.2a) ds^ = -A(r, R)dr^ - 2B{r, R)dRdr - C(r, R)dR^ + R^{de^ + sin^ ^d*^^), 
where 

/ 2vI/\ Y 1 

(2.2b) A= 1-— G^ B = -G-, C = —. 



y R J ' U ' U"^ 

In ( I2.2bl) . the two functions \E'(r, i?) and Y{r,R) are arbitrary (positive) functions 
while 

(2.2c) u' = Y' + ^J-l, 

and the function G is given in terms of a quadrature: 

(2.2d) G(,fi)= r 1 ^(,.)d.. ' 



jji Y{r,a) dr ' Y{r,r)u{r,r) 

In particular, the function \E' represents Misner-Sharp mass. Conditions on Taylor 
developments of the above quantities may be given to characterize complete collapse 
and naked singularity formation, that we summarize in the following theorem. 

Theorem 2.1 ([7|). (1) In the spacetime described by the metric (I2.2al )-( l2.2dl ). 

the singularity forms at the centre r = in a finite amount of comoving time 
if Taylor development of the function Rv? is as follows: 

(2.3) Ru^ = ar^ + ^r'^R + ^rR^ + 5R^ + o{r'^ + R^f/^, (a > 0). 
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(2) Under the above condition, and introduced the Taylor development of the ( reg- 
ular) function G{r, 0) 



(2.4) G{r, 0) = ^r'^-^ + o(r'^-^^ 

then the central singularity is (locally) naked if n = 1, n = 2, or n = 3 and 
^ > a uj cvit where u^nt = 

We also remark that this model possesses nonvanishing anisotropic pressures, which 
are given by Q 

(2.6) Pt = -z KT^iTJTT^-TT^^^+^G 



87TuRG\YdrdR dr OR dR^ 

and in this sense the model can be considered physically more reasonable than Tolman- 
Bondi dust collapsing sphere, where the spacetime is ruled by the pressureless equa- 
tion of state Pr = Pt = 0- 

2.2. Exterior region. In our model, since the internal source is given by (I2.2al ). one 
cannot hope, in general, that radial pressure vanishes at some timelike hypersurface 
S = {(t, r, 6*, 0) : r = rt}. This is a quite restrictive feature of dust cloud collapsing 
model lfT3l or also non vanishing radial pressure models 1 1 1 1, that here arise only as 
very special cases, that happen when ^ = ^(r). Therefore, it cannot be possible 
to consider Schwarzschild vacuum solution as external region, if we want Israel- 
Darmois condition to be satisfied. 

Hence, we perform a junction between the internal source satisfying conditions 
stated in Theorem l2.1[ and the anisotropic generalization of de Sitter spacetime, which 
are a class of spherically symmetric solutions of Einstein equation satisfying the con- 
dition e + Pr = and admitting a particular G4 group of motions (see Q). A coordi- 
nate transformation exists, that brings the line element in the form 

(2.7) ds' = -x{R)dT' + x{R)-'dR' + R'dn^ X(i?) = l-^^, 

thereby obtaining a family of solutions as Misner-Sharp mass M{R) varies. In ||6| it 
is shown that, actually, any spherically symmetric line element (12.11) can be matched 
with a metric of this family at a timelike hypersurface S as before, under the condition 
of continuity of mass only. Therefore, it suffices to choose 

M{R) = m{n,R), Vi?G[0,r,] 

and junction conditions will be certainly satisfied for all (comoving) times t > 0. The 
value of the external mass for bigger values of R depends on the internal spacetime 
at comoving times prior to observation starting, and we will suppose that it is chosen 
such that 

(2.8) lim sup M(i?) < +00. 
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2.3. Energy condition. In order to deal with a physically reasonable class of solu- 
tions, we impose the weak energy condition (w.e.c.) on the energy-momentum tensor 
T of the spacetime. Basically, this means T(t>, v) > for all timelike vectors v, and 
it is easily seen ||3 that it holds in the internal region if, for each r e [0, r;,], the two 
functions of the variable R only 

are nonnegative subsolutions of the same ODE: 

(2.10) 14^™ 

for i? > 0. In particular, for r = r^, the above condition ensures w.e.c in the external 
region also (see ^ eq. (4.6)]). 

Models of internal solutions satisfying ( 12. 3I ). ( 12.81) and w.e.c. can be easily found: 
for instance, the choice 

(2.11) m{r,R)= j -i{s)s^ds+ j x(^)(7^d(T, Y = Y{r), 

Jo Jo 

with 7(s) and xi'^) positive and not increasing in [0, +oo), allows for (12.31) and the 
weak energy condition to be satisfied. As an example, one may take 

in order to satisfy also (12.81) . Note that, with the choice of x = constant, the spacetime 
coincides with the so-called Tolman-Bondi-de Sitter (TBdS), and that's the reason 
why, in |7 1, the models arising from the above choice ( 12.111) are termed anisotropysa- 
tions of TBdS spacetime. 



3. Null geodesics from the singularity 

3.1. Geodesic equations. As stated in Theorem 12. 11 conditions on the metric func- 
tions allow to determine when the central singularity is locally naked. This is made 
by showing the existence of a future pointing radial null geodesic which lies in the 
region R > 2^ and may be traced back to the central singularity. In [21 it is shown 
that, to obtain violations of cosmic censorship in spherical symmetry, one may restrict 
oneself in looking for null geodesic which are radial only. In other words, if a sin- 
gularity is radially censored, it is censored all the way. Nevertheless, if one wants to 
study visibility under a more general point of view, also non radial light rays should 
be taken into account. 

Let K be the affine parameter of the null geodesic. Without loss of generality, we 
will suppose that the geodesic lies in the hypersurface 6 = ^, and then the angular 
components of tangent vector along the geodesic read 



(3.1) 



dK 
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where (. is the (conserved) angular momentum. Let us also introduce a function q such 
that 

dr 1 

(3.2) ^'■■=:r = 

an qRG 

Then, expressing R, 9 and as functions of r instead of k, the equations for null 
geodesies can be given in the following form: 



(3.3a) — = uG 



Y -u\ 1 



tq ' 

u 



(3.3b) 
and 

(3.3c) 



dr 

d^ _ Q d0 _ IqG 

dr ' dr i? ' 

1 dg 1 dG 1 di? C A,. + - 25 , 



2 



g dr G dr i? dr 2G 

_ CA^R-BC,. dR BC,R + C{C,,-2B^n) (^\\ f<f^ ^ 
G2 dr 2G2 \dr ) R 

and A, 5, C are given in (I2.2bl) . Dealing with radial geodesies results in vanishing of 
and so (I3.3al) . which is found imposing that the geodesic is null, becomes an ODE in 
the R function only, which is enough to study at least the existence of corresponding 
pregeodesics. On the other side, when £ 7^ 0, also (I3.3cl) must be studied. 

Before to state and show results, we restrict the analysis hereafter to the case when 
n in ( 12.41) is equal to 3. As it will be cleared in Remark this particular situation 
- that, as seen in the statement of Theorem 12. 11 has a sort of "endstate transition" - 
corresponds to a so-called strong curvature singularity, unlike n = 1,2 cases - at 
least along null radial geodesies. We refer the reader to Remark l34l below and, for 
instance, to |IT1[T71[TH1 for a general insight about strong curvature conditions. 

3.2. Null geodesies existence. Existence of null radial geodesies has already been 
proved in |j|7j, using comparison arguments in ODE. We here complete the analysis 
about the asymptotic behavior of the geodesies finding also an existence result for 
nonradial geodesies. 

Remark 3.1. In the forthcoming Proposition it will be shown that the null geodesic 
equation can be put in the form 

(3.4) ^ = -/(!/) + 

dr r 

where y{r) G E" and / : E" — E, (7 : E""^^ — > E are functions; Introducing 
a new independent variable s such that ^ = — r(s) then, from (13. 4I) . we obtain the 

f^=-f{y{s))-r{s)g{r{s),y{s)), 
fs = Ms), 

whose solutions describe parameterizations of solutions of (13.41) by the parameter s. 
Equilibria of the above system (as s +00) are given by points (yo, 0) where yo is a 
root of f(y) = 0, and searching for these equilibria amounts to search for admissible 
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solutions of the so-called root equation, first introduced in |l3| for the study of radial 
light rays. 

The key point is that these equilibria for the null geodesic equation will be showed 
to be hyperbolic - see e.g. lfT6ll for basic concepts about hyperbolic dynamical systems 
- and that is the reason why the root equation solution's existence is a necessary but 
also sufficient condition for light rays existence. 

Proposition 3.2. Under the hypotheses ofTheorem \2.1\ when , with reference to equa- 
tion (12.41) . n = 3 and ^ > atUcrit, there exists infinite null geodesies emanating from 
the central singularity and escaping from the trapped region {R > 2\E'}. In particular, 
there exists three numbers xi, X2, Xc, depending on a and C,, with Xc > X2 > xi > a, 
such that these geodesies can be divided in the following classes: 

(1) infinite radial and non radial geodesies such that R{r) = xir^ + o(r^), 

(2) infinite non radial null geodesic such that R{r) = X2T^ + o(r^), and 

(3) a radial geodesic (Cauchy horizon) such that Rcij) = XcT^ + o(r^), which 
bounds from above any geodesic in classes ^ and ©, in the sense that 
Rcir) > Rg{r),'ir > 0, for any other geodesic Rg{r). 

Moreover, the nonradial geodesic of case O have finite limr^o+ 4>{''^)' whereas 
nonradial geodesies of case © are such that \imr^o+ (j){r) = — oo. 

Proof Let us consider system (I3.3al )-( l3.3cl) . We are interested in determining exis- 
tence of solutions such that R{r) lies above the apparent horizon Rhir) , which is 
known (see [7| for details) to have the behavior Rhir) = ar^ + o(r^). Then, we first 
introduce a new unknown function z{r) in place of R, which is defined as 



and we will study the system made by the first and the last equation above, in the 
unknown functions {z{r), g(r)) - with g(r) given by (13.21) - since the equation for 
can be decoupled from these. Moreover, note that since i appears in (I3.3al )-( l3.3cl) as 
a factor of the quantity iq only, we will consider iq{r) instead of g(r) as variable. 

We first observe that z{r) and g(r) must be bounded from below by 1 and 0, re- 
spectively, in order for the solution to be physically acceptable. Then, estimates of 
terms involved in (I3.3al) will be performed, under the additional hypothesis that z{r) 
is also bounded from below, and equations (I3.3al) and (I3.3cl) become, with the further 
position 



(3.6) 





(3.7b) 



(3.7a) 



dz 

dr 
^dq 
dr 



r 



1 




1 



w{z,iq) +gi{r, z,iq), 
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where 




(3.7c) 


w{z, £q) 


(3.7d) 


k{z) = - 


(3.7e) 


h{z) = - 


(3.7f) 


k{z) = - 


and gi{r,z, 





2z^ + ^/iTWl^z^ -UOZ+ ^/iTW)^^ 
uo^ - IQz^uo - 2z^ 



2z\iU + z3) 
u + z^ 



5 

z 

z^ + Au 



z{uj + z^) 

1, 2 are continuous functions. 
Let us check equilibria of (I3.7al) - (l3.7bl) . using the idea from Remark ISTTl The 
following cases happen. 

(a) £q = 0, z such that w{z, 0) = 0; 



First, in both cases it can be seen that ((71 (r, z, iq),iqg2{r, z, iq)) is regular in a 
neighborhood of (0, z, £q) where {z, iq) is an equilibrium as above. Let us now dis- 
cuss the character of these equilibria. 

In case (a), the situation goes as follows. With reference to the notation of Remark 
13.11 the Jacobian computed in the equilibrium (z, 0) is given by 



(3.8) J(z,o)f 



k{z) 
and one finds that there exists a root function 

2z^ + z^ 



(3.9) piz) 
such that 



z-1 



(3.10a) w{z, 0) = ^ p{z) = uj, 

(3.10b) ^(.,0) = -^p'(.). 

Moreover, k{z) > if and only if > p{z) := (5 + v^27)-2;^. The curves p(z) 
and p{z) are depicted in Figure [TJ For the acceptable range of u, that is > c^crit 
(see Theorem 12.1 1) , there exists a unique value Zi that satisfies (I3.10al) and such that 
k{zo) > 0. From (IB.lObll . we have that also ^(2:1,0) > 0. This means that the 
Jacobian ( 13.81) has two positive eigenvalues and therefore we have part ([T]) of the 
claim, with an infinite number of radial geodesies tending to this equilibrium point, 
for every £ > 0. Note that this situation (HI) cover also radial geodesic existence, that 
is when i = 0. Using (13.61) . the curves R(r) behaves like R{r) = azfr^ + o{r^). 

If UJ > uJcTit, from Figure [T] we see that there exists another values of z = Zc 
satisfying p{zc) = uj- In this case we have two negative eigenvalues of (13.81) . With 
reference to the system put in the form (13.51) . the equilibrium is again hyperbolic, but 
the stable manifold has dimension 1, and therefore a single null geodesic tending to 
the equilibrium exists. Actually, it can be shown that this geodesic is radial, since 
it can be seen that at least a radial null geodesic tending to the equilibrium exists. 
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Zl 



Zc 



Figure 1. The curves p{z) ( I3.9I) . a{z) (I3.12I) . and fi{z) intersect each 
other at the local minimum of p{z), attained for z = ^-y^, p{z) = 
cUcrit (see Theorem 12 .11 1. 

Indeed, radial geodesies are solutions of the ODE ( I3.7al) under the assumption i = 0. 
We can easily see that z = Zcis an equilibrium of (I3.7al) . such that there exist a unique 
solution tending to this equilibrium. Since radial geodesies are solutions of (I3.7al) - 
(I3.7bl) with dq = 0, we conclude that the (unique) solution of (I3.7al) - (l3.7bl) must be 
radial. The corresponding x = azi represents Cauchy horizon direction. 

We also note that the corresponding function Rc{r) must bound from above any 
other geodesic Rg{r). Of course, this must be true for any other radial geodesic near 
r = 0, and therefore for any other value of r by uniqueness of Cauchy problem solu- 
tion. Let us show that this also holds true for nonradial geodesic. By contradiction, 
let Rg{r) be a nonradial geodesic, such that its projection on the (r, R) plane bounds 
from above Rc{r) for some r^, that is Rg{r^) > Rc{r*). Since Rg{r) is a subsolution 
of radial null geodesic equation, and Rc{r) is a solution of the same ODE, therefore 
we can find a radial null geodesic Rn{r) with initial data Rn{r^) g]-Rc(^*), 
and trace it back until _R„(0) = 0, obtaining a radial geodesic _R„(r) > Rc{r) inaright 
neighborhood of r = 0, which as said before is absurd, and claim Q is complete. 

Let us now discuss case (b) to get claim Q of the Proposition. 
We first observe that this time it must be k{z) > otherwise we have no equilibria. 
Arguing as in the previous case, we find that the determinant of the Jacobian in the 




(3.11) 




The root function in this case is given by 
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which has the properties 



(3.13a) w{z2, = ^ ^(^2) = ^, 



where 7 is some non vanishing constant. The situation is depicted in Figure [T] For 
any uj > Ua-it there exists a unique value of Z2 satisfying (I3.13al ) with A; (-22) > 0. 
Using ( I3.13bl ). we get that the quantity O.l II) is negative and then there exists one 
and only one positive eigenvalue of the Jacobian, which results in the existence of 
infinite non radial null geodesic such that the corresponding R{r) behaves like R(r) = 
az^r^ + o(r^). Indeed, the stable manifold of the system put in the form (13.51) has 
dimension 2. 

Notice that, in this case, there are no other choice available since the other root of 
the equation a{z) = uj lies in the {A; < 0} region. 

To complete the proof, let us now analyze the behavior of the angular function 0(r). 
Using (I3.3bl) . and the above estimates, 

(3.14) ^<P^JKz)\q 



dr \ z J r 

For a non radial geodesic from case ([T]), g ^ as r — 0+, and then, from (I3.7bl ). 
^ = k{zi)^, from which q{r) = r''(^^\ and then, using (13.141) . 0(r) has a finite limit 

as r — > 0^. 

For a non radial geodesic from case (13) conversely, q tends to a finite positive value 
and then the function 0(r) goes like kr^^, which determines a negative diverging 
behavior of the function 0(r) . □ 

Remark 3.3. These results are completely consistent with the analysis carried out in 
lfT2l fT4l. which have established a double topological nature of the naked singularity 
in the particular case of dust solution under self-similarity assumption. In particular, 
case (UJ) corresponds to a region of the singularity foliated by a 2-sphere, whereas 
case (13) corresponds to a topologically pointwise singular region. 

Remark 3.4. If we look at the behavior of the quantity k'^Ra/sK'^K^ along null radial 
geodesies, where Ra/3 is Ricci tensor, K"" = ^ is the tangent vector of the geodesic 
with parameter k, we get 

(3.15) eR.„K'K^^^-^(K«r-. 

In view of the above theorem, equations ( 13.21 ) and ( I3.3al) ensure that there exists 
limfc^o+ then tends to a finite nonzero limit value. Since both G and \1/ 
goes like r^, and Y, u tends to a finite nonzero limit along the geodesies, the quantity 
in (13.151 ) goes to a finite nonzero value. It follows that these null geodesies terminate 
in a strong curvature singularity in the sense of Tipler ifTTI . 
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4. Faraway visibility of the singularity 

4. 1 . Size of the singularity. In the following, we will study the behavior at infinity of 
null geodesies emanating from the singularity. We have already seen that there exists 
an infinite number of light rays escaping from the trapped region; These geodesies 
arrives at the boundary of the collapsing sphere of anisotropic matter with some value 
R{rh), and they have to be continued by studying the behavior of null geodesies in the 
external region (12.71) escaping from the boundary and such that R = R{rh). 

The equation for null geodesies - that, without loss of generality, we will suppose 
to lie in the subspace {9 = f } - is given by 



where k denotes as usual the affine parameter and ^ and uj are constant of motion. 

As explained in fTT|, since the metric (12.71 ) does not depend on T we can uniquely 
associate to an observer R = R^aa orthonormal basis of the tangent space to a point 
(T, Rq, 6, (f)) of the outer region (i.e. such that Rq > 2M{Rq)), which will be denoted 
by 

{e(T) = x{Ro)~^^^dT, e^R) = xiRoY'^dn, e(9) = R^^de, e(<^) = {R^smdy^d^] . 

Then, if we consider a null geodesic with affine parameter k emanating from the 
singularity passing for R = Rq, and denote by = its tangent vector, this 
observer measures an angle 6 between the light ray and the radial direction equal to 



(4.2) S = arctan 



df^c^e^R)!^^ u\j Rl- {i/uyx{Ro) ' 



The above quantity depends on the geodesic, and therefore the supremum made among 
the set & of all singular geodesies detected at R = Rq gives a measure of the singu- 
larity detected by the observer. One can conceive the righthand side above, fixing 
-Ro' as a function in (i/to) which results to be increasing. Therefore, the "size" of the 
singularity detected by the faraway observer is related to the quantity 

:= sup — , 

e ^ 

which can be regarded fV2\ as a sort of impact parameter. 

Using this definition, it is not hard to prove an extension of the result already proved 
for the particular case studied in lfT2l . 

Proposition 4.1. If the boundary of the interior region r = Vh is sufficiently small, 
then the non radial geodesies of largest impact parameters emanating from the singu- 
larity are such that the angular function — * — oo m the approach to the singularity 
(case (121) of Proposition 13. 21 see also Remark \T3\l . 

Proof. Of course, we will consider only nonradial geodesies, when £ ^ 0. 
Using continuity of the metric along the junction hypersurface, one gets 

dr 



dT dTdi? dTdr 



Y dR ^ 

TT + ^ 

Xu ar 



dfi;' 
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and using (14.11) . together with (I3.2l) - (l3.3bl) one finds the angular frequency uj and 
therefore the impact parameter 6, which is given by 



(4.4) h = sup 

6 



R{r) 



Y{r, R{r))^u\r,R{r)) + (£g(r))2 - u^{r, R{r)) 

With reference to the cases listed in Proposition l3.2l let _Ri (r) be a nonradial geodesic 
from case (O, and let R2{r) a (necessarily non radial) geodesic from case @. It is 
easy to verify that, as r ^ 0"*", the quantity in round brackets in (14.41) tends to along 
i?i(r) and tends to a finite nonzero value along R2{r). Then, if rb is sufficiently small, 
we can suppose that this quantity, along R2, bounds from above the same quantity 
computed along Ri until r = Vb- □ 

4.2. Redshift. Following the frequency shift z between a source and an observer 
respectively located at events Pi and P2, connected by a null geodesic with tangent 
vector ~K w.r.t. the affine parameter, is defined as 

Op, ("/?, U(s)) 

(4.5) 1 + 



gp^{ K ,M(o))' 

where u^s) and are the 4-velocities of the source and the observer respectively. 

To determine the redshift associated to the singular geodesies, the numerator above 
will be replaced by a limit expression as Pi approaches the central singularity. Let us 
therefore consider a singular geodesic passing at event P2 in the exterior space such 
that R{P2) = Rq- Using the basis (14.31) . together with (14.11) . it is easy to see that 
the denominator in ( 14.51) is given by — ^(-Ro)^^^^^ which under the hypothesis ( 12.81) 
is always finite also if Rq tends to infinity. On the other side, the velocity field in the 
interior spacetime is given by —u{r, R)dR, and using (I3.2l) - (l3.3al) one finds that the 
absolute value of the numerator in (14.51) is given by 

1/2 



(4.6) lim 



0+ R{r) 



u{r, R{r)) 
q{r) 



2 



+00, 



which results in an infinite frequency shift. Actually we can observe that calculation 
of (14.61) above seems a feature of all singular spherically symmetric models, at least for 
nonradial geodesies, when i ^ 0. In the case under our study, moreover, (14.61) holds 
also for radial geodesic since, using Proposition ^ .2l and equation (I2.2cl ). u{r, R{r)) = 
Uq 7^ 0, and both R(r) and q(r) are infinitesimal. It appears also to be possible the 
study of the limit in (14.61) in a even more general framework, and work in this direction 
is in progress. 

All in all, these models shows the same feature described in ||4| for dust solutions: 
in the case corresponding to a strong naked singularity (along radial geodesic^), null 
geodesies are infinitely redshifted. Therefore, it can be said that there exists a form of 
weak censorship for such class of solutions, at least at the classical level. Of course 
however, quantum phenomena might change this scenario, as suggested by some au- 
thors (see irnii.O). 



^As pointed out in |4|, it can be possible that the strength of the singularity may have a directional 
behavior, that is the singularity is weak - in the sense explained in Remark l34l - along some geodesies, 
and strong among others. 
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5. Discussion and conclusions 

In this paper we have studied visibility of the naked singularity for a faraway ob- 
server, investigating global behavior of null (radial and nonradial) geodesies in a class 
of spherically symmetric singular spacetimes. The interior part belongs to a wide class 
of anisotropic solutions; as known from the analysis carried out in [7 1, the endstate of 
the central singularity is determined by the lower order term n of Taylor development 
of a quantity determined by the metric. In this paper we considered the case n = 3, 
which has a special interest since (i) the endstate is determined by the value of a cer- 
tain parameter ^/a, that causes a sort of "phase transition" between naked singularity 
and black hole situations, and (ii) radial geodesies terminates into a strong curvature 
singularity, as observed in Remark 

The exterior part belongs to the class of anisotropic generalization of de Sitter solu- 
tion discussed in [5 1. The form of the metric closely resembles the Schwarzschild line 
element, but the mass here depends on the areal coordinate R since a constant mass 
would not properly match an interior solution with non vanishing radial pressures. 
Indeed, if Misner-Sharp masses for both solutions (interior and exterior) coincide at a 
hypersurface r = const., this ensures Darmois junction conditions to hold along this 
hypersurface. 

Existence of nonradial geodesies is crucial to investigate global visibility. The sys- 
tem of ODE satisfied by the nonradial geodesies does not decouple and the most 
suitable way seems to study equilibrium points of this system. It turns out that there 
are two values xi < X2, depending on the data of the spacetime, such that there is an 
infinite number of nonradial geodesic R{r) emanating with "direction" xi (i.e. such 
that lim^^o = xi) and infinite nonradial geodesies with direction X2. This 
result is not in contrast with previous results regarding self similar collapse [|14I where 
one has infinite geodesic with direction xi vs. only one geodesic with direction X2. In- 
deed, as the proof of Proposition l3.2l points out, the stable manifold of the equilibrium 
has dimension 3 in the first case and 2 in the second one. Self-similarity assumption 
simplifies the geometry of these stable manifolds, shrinking them by one dimension. 

It must be remarked that a crucial fact is that ODE system can be brought in the 
form (13 .4I) . where the equilibria are all hyperbolic. This is a feature of all collapsing 
matter models known in literature, which excluded pathological situations such as the 
ones outlined in [13|, and explains why the so-called root equation method repeat- 
edly exploited in previous works - based on the existence of a certain limit which in 
principle cannot be taken for granted (see for instance |I21lll[l2J) - leads all the same 
to the conclusions one would expect. 

The analysis of the ODE system allowed us to determine a sort of "impact parame- 
ter" that characterizes the "size" of the singularity as it is seen by a faraway observer. 
It turns out that photons emanated by the naked singularity are infinitely redshifted 
also when pressures are present, as has been already shown (see [4|) for the - pres- 
sureless - Tolman-Bondi collapse. 
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